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Topological superconductors supporting Majorana Fermions with non-abelian statistics are 
presently a subject of intense theoretical and experimental effort. It has been proposed that the 
observation of a half-frequency or a fractional Josephson effect is a more reliable test for topo- 
logical superconductivity than the search for end zero modes. Low-energy end modes can occur 
accidentally due to impurities. In fact, the fractional Josephson effect has been observed for the 
semiconductor nanowire system. Here we consider the ac Josephson effect in a conventional s-wave 
superconductor- normal metal-superconductor junction at a finite voltage. Using a Floquet-Keldysh 
treatment of the finite voltage junction, we show that the power dissipated from the junction, which 
measures the ac Josephson effect, can show a peak at half (or even incommensurate fractions) of 
the Josephson frequency. The ac fractional Josephson peak can also be understood simply in terms 
of Landau-Zener processes associated with the Andreev bound state spectrum of the junction. 



Introduction Topological superconductors pQ are 
promising candidates for the practical solid state real- 
ization of Majorana Fermions (MF) 51 [TT]. The MFs 
are predicted to occur as zero-energy bound states at- 
tached to defects and have received much attention re- 
cently p~2j [13] , both due to their predicted non- Abelian 
statistics [T4UT7] and their potential application in topo- 
logical quantum computation (TQC) [T3l [T8H2T] . A sim- 
ple topological superconducting (TS) system supporting 
MFs, which has attracted considerable experimental at- 
tention [13] , consists of a spin-orbit coupled semiconduct- 
ing system in a magnetic field placed in contact with an 
ordinary superconductor [7HTT] . MFs at the ends of such 
a wire have been predicted to produce a zero-bias con- 
ductance peak [3] [TTJ [22]. In fact, recent experiments 
[23H25] measuring the tunneling conductance in the semi- 
conducting wire system suggest the existence of MFs. 

While the zero-bias conductance peak observations 
[23H25] are encouraging, peaks in the tunneling spectrum 
arising from impurities are difficult to completely rule 
out. The end MFs appear in the TS system as a result of 
a Pfafflan topological invariant associated with the bulk 
[TJ [4] [26] of the wire. In a TS system in the ring geom- 
etry, with a weak-link (shown in Fig. [TJa)) connecting 
the ends, the bulk topological invariant leads to topo- 
logically protected crossings of localized Andreev bound 
states (ABS) as a function of the flux penetrating the 
ring [4] [9]. The crossings of the ABS of an SNS junc- 
tion can, in principle, be inferred using the ac Josephson 
effect with a finite DC voltage V across the junction so 
that the phase varies in time as (j>(t) = 2tt^ = Qjt. Here 
Qj = 2eV is the Josephson frequency, where we have set 
h = 1. If the applied voltage V is small enough, the time- 
evolution of the junction can be considered to be adia- 
batically following the ground state of the SNS junction 
with a phase difference <fi. The current in the junction 
I(t) = V- 1 ^ = ^ (where E is the total energy) then 
varies in time corresponding to the ABS energies. If a 
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FIG. 1: (Color online) (a) Low-density normal region with 
chemical potential \i in between two s-wave superconducting 
(SC) leads with phase difference cj). The normal region is long 
enough to support confined ABSs with energy less than the 
SC gap Ao in the SC leads. (b)ABS energy (E) spectrum 
in the normal region as a function of cj) for A = 0.4 meV 
and fi = 0.13 meV. The ABSs in the conventional supercon- 
ducting junction (solid and dash-dotted lines) show a weak 
avoided crossing at <j> = 7r[321 133] . while the ABSs in the 
topological superconducting states (dashed and dotted lines) 
show a protected crossing [4] [9] 1101 128] . 

pair of ABSs (shown as dashed lines in Fig. []Jb)) cross 
zero energy, evolution of the state following the ABSs 
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as <p varies from to 2tt leads to a state with a differ- 
ent energy at the end of each period t = r = 2tt/VLj. 
This results in a current-phase relation with a compo- 
nent I ((f)) ~ sin ^, which is at double the period of the 
conventional Josephson effect [4j [9j [10, 27-29 , and is 
referred to as the fractional Josephson effect. Here (f> 
is the gauge invariant phase across the superconducting 
tunnel junction and is given by (f> = 27r<l>/ ( I>o, where 
is the magnetic flux through the loop and $o = hc/2e. 
However, at finite temperature or in the presence of a fi- 
nite density of fermionic states, the excited ABS state for 
(f> > 7r can relax to the lower energy state so that frac- 
tional Josephson effect in TS systems can be observed 
only above a finite frequency [27 . For conventional su- 
perconductors without such zero-energy crossings of the 
ABSs (solid lines in Fig. [I](b)), the current in the adia- 
batic and low-temperature limit is 27r-periodic. 

The periodically varying current /(£), for a junction at 
constant voltage V, will lead to radiation with a power 
spectrum P(oS) oc< |/(cj)| 2 > with peaks corresponding 
to the Fourier components of the current-phase relation 
I ((f)). TS systems in such configurations have been pre- 
dicted to show peaks at half the Josephson frequency 
i.e. at uj — flj/2. The power-spectrum radiated from 
junctions is somewhat difficult to measure and it is more 
convenient to look for absorption of microwaves by an 
SNS junction leading voltage steps called Shapiro steps. 
The Shapiro steps occur at voltages where the Joseph- 
son frequencies of the corresponding voltage matches the 
applied frequency [30]. Indeed, the fractional current- 
phase relationship has been observed in the semiconduc- 
tor nanowire system in the form of double voltage Shapiro 
steps, strengthening the evidence for the topological su- 
perconductivity in these systems [31] . 

In this paper, we show that, even non-topologial super- 
conducting systems, such as a conventional SNS junction 
constructed out of s-wave superconductors can produce 
an ac fractional Josephson (or Shapiro) effect in P(cj), 
if the applied voltage V happens to be larger than an 
avoided crossing gap (shown in Fig. [I]) that can acciden- 
tally become small, even though the current-phase rela- 
tion of such a system would be 2tt periodic. Thus, similar 
to zero-bias tunneling signature for MFs, the fractional 
Josephson effect can also occur accidentally in conven- 
tional systems. For our study, we will ignore the effects of 
Coulomb interactions. In the high transparency param- 
eter regime that we study one can expect these effects to 
be small and only effectively renormalize the parameters 
such as position of the crossing of ABSs as a function of 
chemical potential. The conventional SNS junction only 
serves as a model to study the fractional Josephson ef- 
fect. Similar physics is expected to hold in other systems 
where pairs of ABSs can cross zero-energy as a result of 
tuning of parameters. 

To see how an SNS junction can give rise to a frac- 
tional ac Josephson effect, we consider the Hamilto- 



nian of an s- wave SNS junction shown in Fig. [TJa) un- 
der a finite voltage V, which is written as h((f>(t)) = 
f dx^(x)h[(j)(t)]^(x) where ^f(x) = (^(x), ifj\(x)) is 
the fermion operator in Nambu spinor notation and 
(f)(t) = Qjt is the time-dependent phase difference be- 
tween the left and the right superconducting contacts. 
Here h is the Bogoliubov-de Gennes Hamiltonian, given 

by 

h[<f>] = [(-V 2 -n)T z +A {(e(x+W)+e(W-x)e i<i >)T + +h.c}}, 

(1) 

where t z ^+- are the corresponding particle-hole Pauli 
matrices. In Eq. [I] Ao is the superconducting gap in 
the S region, 2W is the width of the junction, \i is the 
chemical potential, which controls the density of elec- 
trons in the normal region and Q(x) — 1 for x > and 
Q(x) = for x < 0. Such SNS junctions have been 
shown to have ABS spectra with avoided crossings con- 
trolled by the transparency of the barrier [32]. It has 
also been shown that Landau-Zener (LZ) crossing pro- 
cesses can play an important role in determining the dc 
current response [33 in conventional SNS and the finite- 
frequency response [27] in TS junctions. In the case of 
a high transparency interface with reduced Fermi-energy 
in the N region leads to an ABS spectrum shown in Fig.[l] 
(solid lines) with a gap near (f> = 7r, which is controlled 
by the chemical potential /i. The phase-dependent super- 
current operator I ((f)) is obtained as a derivative 

7(4>) = cV#]. (2) 

For conventional SNS junctions, one expects the ABS 
occupation to stay in the ground state for voltages 2eV <C 
Eg, where E g is the minimal gap in the junction. How- 
ever, typical high-transparency and near resonant SNS 
junctions show only weakly avoided crossings (as seen in 
Fig. [T]) , so that so that 2 eV can be equal or larger than 
Eg leading to a violation of the adiabaticity requirement. 
In this case, the state of the system can crossover from 
one state to the other by LZ tunneling leading to a sit- 
uation analogous to the ABS in the TS case (shown by 
dotted lines in Fig. [I]), resulting in a fractional ac Joseph- 
son effect even in the conventional SNS junction. For 
simplicity, we consider values of \i and W, where the gap 
in the spectrum (between the lowest positive energy ex- 
citations) near (f> ~ is much larger than the energy of 
the lowest energy excitation near <f> = tt (as in solid lines 
in Fig. [TJb)). In this parameter regime, a range of volt- 
ages V allow the dynamics of the lowest energy states to 
remain confined to the lowest energy pair of states. 

The power-spectrum P(uj) dissipated by the fluctuat- 
ing current I(t) in the Josephson junctions with a weak 
applied DC voltage V, can be calculated [35 using the 
expression 

P(w) = P c (oj) + Y,\In\ 2 S(ojT-2wn), (3) 
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where P c (uo) is the connected part of the power spectrum 
defined as 
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(i(h)i(t 2 )) c , (4) 



with (/(^)/(i 2 )) c = </(t 1 )/(t 2 )) c - 

fourier transform I n is defined as 
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Here we have generalized the expression appropriately to 
include discrete time-translation invariance of the Hamil- 
tonian (i.e. t — » t + r). The time-periodic expectation 
value of the current in Eq. [5] gives rise to singularities in 
the power-spectrum in Eq.[3| 

For simplicity, let us first consider the case where the 
current operator I(t) is strongly coupled to ohmic dis- 
sipation from a resistive shunt so that quantum fluctu- 
ations of the current I(t) are suppressed. The current 
operator I(t) can then be replaced by a classical random 
variable I(t) —> I(t). In this classical limit, the current 
I(t) depends on the quasiparticle occupation of the SNS 
junction according to 



I(t) = s(t)i (t), 



where zq(£) is the current of the ABS in the lower energy 
state and s(t) is a random variable such that s(t) = — 1 
in the higher energy state in Fig. [I] while s(t) = 1 in 
the lower-energy state. For now, we have assumed that 
the fermion parity in the junction is fixed, so that the 
junction has only two states. The state variable s(t) is 
assumed to be constant, except near the avoided cross- 
ings where <j)(t) = 2eVt crosses (2n + l)7r, where an LZ 
crossing can transfer s(t) between the values ±1 with 
some probability. The evolution of s(t) can then be de- 
scribed by a classical Markov process with a transition 
probability matrix P [s(t = nr) — » s(t = (n + l)r)]. Us- 
ing the transition probabilities, the connected part of the 
power spectrum for the current fluctuations according to 
Eq. [4] is found to be 
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where io(i) = J2 m i m ^ l2rneVt and a, b are the conditional 
transition probabilities at each LZ crossing at time t = nr 
so that a = P [s(t = (n + l)r) -l\s(t nr) +1] and 
b = P [s(t = (n + l)r) = +l|s(t = nr) = -1]. At high 
voltages V, which are much higher than the avoided 
crossing, where the LZ crossing probabilities a ~ b ~ 1, 
it follows from Eq. [7] that the radiated power P(uS) has 
a peak(from the vanishing denominator) at uo ~ 7r/r = 



flj/2 that is characteristic of the fractional ac Josephson 
effect. In addition, for asymmetric transition probabili- 
ties a ^ b the power-spectrum in Eq. [3] contains a train 
of singularities at uo = 2ir/n with strength 
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For large voltages V, one expects the state of the 
junction to become highly excited and possibly change 
fermion number while remaining quantum coherent over 
a few periods r. Therefore, we consider an SNS junc- 
tion described by the time-dependent BCS Hamiltonian 
h[(j)(tj\ (Eq. [l]), while being weakly tunnel-coupled to 
a fermionic bath. The Keldysh time-contour Green- 
function [34H36] for the system including the time- 
dependent Hamiltonian h[<p(t)] is written in the form of 
a Dyson equation 



G = g + gh[<t>(t)]G, 



(9) 



where g is the Keldysh time-contour Green-function for 
the system excluding the time-dependent part of the 
Hamiltonian(i.e. h[cj)(t)]). The Green function g is thus 
an equilibrium Green function is written as 



( 6 ) 9r,a{u) 
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uJ±iT , 



±2iT 

o 2 + r 2 



n F (±co), (10) 



where T is the imaginary part of the self-energy (assumed 
to be independent of uj) arising from coupling to the 
fermionic bath and np(uj) = 1+e ("L/ T ) is the fermion dis- 
tribution function |35l |36| • Using the Green function in 



Eq.[l0]in the Dyson equation Eq.[9]one finds the retarded 
Green function to be 



G RA (t,t') = Ti#(±(t-t'))U(t,t'), 



(11) 



where 0(t) = e - rt G{t), and U(t,t') = Te^ti ^(tO. 
Expanding out the Dyson equation Eq. [9] [55] . and us- 
ing the relations U{t\,t2)h{t2) = —idt 2 U{t\,t2) together 
with h(t\)U(ti,t2) = idt 1 U{t\,t2), the distribution func- 
tions are found to be 
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T2iG R h F {±{t-t'))G / 
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where npit) is the Fourier transform of np(uj) In the 
T ^> V limit, we can approximate hp ~ S(t — t f ) and 
write 



G < ' > (t,t / ) = ±2ie- v \ t - t '\U{t,t'). 



(13) 



Using the Floquet theorem [37], the unitary time- 
evolution operator £7 (£1,^2) corresponding to the time- 
periodic Hamiltonian h[<p(t)] can be constructed in terms 
of Floquet states ip\(i). The Floquet states ip\(t) are 
defined at t = to be eigenstates of U (r, 0) using the 
relation J7(r, 0)(f\(0) = e l€xT (p\(0). Here 6a are the 
Floquet energies corresponding to <p\(t). The definition 



of (p\(t) is then extended to all time using the relation 
ip\(t + r) = ip\(i) = e- iext U(t,0)(p x (0). The unitary op- 
erator [/(£, 0), which is required for calculating the Flo- 
quet states, can be computed numerically by solving the 
time-dependent Schrodinger equation corresponding to 
Eq. [I] The unitary time-evolution operator can then be 
conveniently written in the basis of the Floquet states as 

U(h,t 2 ) =J2e- i€x(tl - t2) Vx(tM(t2). (14) 



Restricting to the two ABS state case, the Floquet 
energies, e^, are determined from the eigenvalues e zexT 
of the unitary matrix U(r, 0). Since the underlying BCS 
Hamiltonian in Eq. [I] is particle-hole symmetric, U (r, 0) 
is restricted to the form 



so that 



U(t,0) 



2eV 



cos Ve° 
sin 0e 1, 



-sm6e- i(f)1 
cos 6e-^° 



(15) 



e\ = ±— — sgn(cos 0o ) cos 1 {cos^ cos 0o}- (16) 

Z7T 

The connected correlator (I(ti)I(t2)) c defined in the 
context of Eq. [4j which turns out to be a combination 
of the retarded and advanced response functions and the 
shot-noise correlator [34, 36 for the current /(£), is writ- 
ten as 

{i{h)i{t 2 )) c = Tr[I(t 1 )G > (t 1 ,t 2 )I(t 2 )G < (t 2 ,t 1 )}, (17) 

where I[(j)(t)] = d<f>h[<fi\. The connected part of the power- 
spectrum of current fluctuations can be computed using 
Eq.[l7|in Eq.[4| which simplifies after substituting Eq.[l3j 
to 



2r|J AlAa (n)| 2 

r 2 + (c^-cD n ) 2 



(18) 



where uj n = uj - e Al + e\ 2 + nVtj and IaiA 2 (X) = 
dte inQjt ((p Xl (t)\I(t)\(px 2 (t))- It is clear from the form 
of this equation that in the limit 2eV ^> T, one can 
expect the power-spectrum P c (uj) to have peaks at the 
differences in quasi-energy uj ~ e\ 1 — e\ 2 . In the LZ limit 
where the voltage V is much greater than the avoided 
crossing frequency so that — >• tt/2 in Eq. 16 



check from Eq. [16] the difference in eigenvalues 



ei — eo 



2eV 
~2~ 



one can 



(19) 



giving rise to the fractional Josephson peaks in the 
current-noise spectrum. This argument can be verified 
(as shown in Fig.[2| by numerically evaluating the power- 
spectrum in Eq. [l^] corresponding to the ABS spectrum 
in Fig. [I] Interestingly, we find that in the intermediate 



,3, 

c 



.E/V = 0.5 

g dc 

E N A = 0.8 

g dc 

. E /V. =2 

g dc 

.E/V, =4 

g dc 



-E/V. 

g dc 


= 0.5 




... E /V. 

g dc 


= 0.8 




— E/V. 

g dc 


= 2 




— E /V. 

g dc 


= 4 





2 2.5 
V„/(o 



3.5 




0.25 



0.5 



o/Q T 



0.75 




1.25 



FIG. 2: Power dissipated P c (cj) as a function of frequency 
uj/Qj for different values of the ratio of the Andreev gap E g 
(i.e. half the avoided crossing in the ABS near <fi = tt) to the 
applied voltage Vdc> The ratio E g /Vd c is varied by changing 
li. The power-spectrum P c , for small E g /Vd c , shows a peak 
at uj — Qj/2. In contrst, P c {oj) shows only the conventional 
Josephson peak at uj = Qj for larger values of E g /Vd c - Here 
P c (oj) is only the connected part of the radiated power spec- 
trum (see Eq. |3j). The total power P(uj) may contain addi- 
tional contributions at integer frequencies uj ~ nQj. The inset 
shows the Shapiro kink response f(Vd c ) = (hotRi — Vd c )/Ri 
of the DC voltage Vdc as the bias current I to t changes. The 
function / shows peaks at both the conventional voltage 
Vdc = huj/2e and the doubled Shapiro voltage Vdc = 2huj/2e 
as a function of the ratio Vdc/oj. Here Vdc — 1 /xeV and fi 
ranges from 0.12 meV to 0.15 meV. 



regime between the fractional and conventional Joseph- 
son effect, the power-spectrum P(oS) can have peaks at in- 
commensurate frequencies. Additionally, the terms pro- 
portional to I n in Eq. [3] together with other ABSs in the 
conventional as well as TS junctions, which follow the 
conventional Josephson phase periodicity are expected to 
lead to the peaks at the conventional Josephson frequen- 
cies. However, it is possible for the fractional Josephson 
contribution to dominate over the conventional one as 
appears to be the case in recent experiments [3T] . 

Shapiro steps: The time-dependent current in a finite 
voltage biased SNS junctions can be detected more easily 
through the detection of Shapiro voltage steps across the 
Josephson junction [30]. As discussed in the appendix, 
the Shapiro step experiment consists of injecting ac power 
into a resistively and capactively shunted SNS junction 
through the application of an ac voltage and measuring 
the response in the DC voltage across the SNS junction, 
which is biased to be in the over-damped regime (as in 
the ac fractional Josephson effect) so that it always has a 
finite DC voltage Vdc 7^ across it. Therefore the Shapiro 
step experiment is the reciprocal measurement of the ac 
Josephson effect in essentially the same set-up (which is 
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shown in Fig. [3] of the appendix). In a conventional (low 
transparency) SNS junction where Vdc is much smaller 
than the minimum gap in the junction, the application 
of an ac voltage with frequency uj leads to a plateau in the 
dc voltage across the SNS junction as a function of I tot 
at Vdc — oj/2e. The range of I tot over which this plateau 
exists is found to be proportional to the amplitude of the 
applied ac voltage v aCy o (as reviewed in the appendix). 
For a finite temperature or otherwise noisy classical SNS 
junction, such as the one described by Eq. [6] earlier in 
the text, we expect that the SNS junction does not carry 
long time classical correlations. Therefore the response 
to v ac ^o should be analytic in the v aCy o — ^ limit and also 
be time-translation invariant. As a result the response 
of Vdc scales as v 2 ac , which is different from the plateau. 
To compare the small v ac ,o response with the plateau it 
is convenient to define the function 



f(V dc = I tot R 1 ) = 



ItotRl — Vdc 

Ri ' 



(20) 



which in the case of a plateau shows a kink near specific 
values of Vdc that resembles the derivative of the peaks in 
the power-spectrum of the ideal SNS junction in the limit 
of small uj (see appendix for explanation). In the classical 
fractional junction case, where the power- spectrum has 
a peak at uj — eVd c (see Eq. [7]), we find a similar plateau 
with width in current I tot proportional to v ac ,o, when- 
ever v aCy o is large compared to the noise induced broad- 
ening of the power-spectrum peak. In the small v ac ^ 
limit we argue that this plateau broadens out into a kink 
(henceforth referred to as a Shapiro kink ) whose weight 
is proportional to v 2 ac and whose profile resembles the 
derivative of the broadened power-spectrum peak. The 
resulting shape for f(Vd c )/Vac,o 1S related to the response 
of the noise correlator of the current to shifts of the phase 



f(V di 



u ac,0 



dtdt'dt F(t-t ,t-t') 



6(I sc (t)I sc (t )) 
Hd{t') ' 



(21) 



which depends only on the linear response properties of 
the junction and F{t — t$,t — t f ) written out explicitly in 



the appendix (see Eq. 49 ) is a response kernel of the ex- 



ternal RC circuit. The resulting kink profile function / 
calculated using Eq. [21] as seen from the inset in Fig. |2j 
shows a kink corresponding to the derivative of the dissi- 
pated power-spectrum from the ac fractional Josephson 
effect. 

Conclusion: We have shown that high transparency 
conventional SNS junctions with weakly avoided zero- 
crossings [32] can show a peak at uj ~ Qj/2 in the power- 
spectrum thus showing a fractional ac Josephson effect 
even at voltages V <C A. In the case of coherent dynam- 
ics of the SNS junction, we find peaks in the radiation at 
a frequency between flj/2 and Slj, which we conjecture 
could lead to Shapiro kinks at incommensurate values. 
The fractional Josepshon effect arises here as a result of 



LZ processes resulting from a break-down of adiabatic- 
ity near the avoided crossing and exists for voltages that 
are larger than the gap E g in the Andreev spectrum. In 
contrast the fractional Josephson effect in TS systems 
arises from protected zero-energy crossings in the ABS 
spectrum [4j [28] and therefore would theoretically exist 
up to V —> 0. In practice, damping from fermionic baths 
restricts the observation of the fractional Josephson ef- 
fect to finite (and often large) voltages V. Therefore, 
the ac fractional Josephson effect at finite voltages, sim- 
ilar to zero-bias conductance peaks can arise both arise 
in conventional superconducting systems as a result of 
accidental fine tuning. 
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FIG. 3: Schematic circuit diagram for studying the ac Joseph- 
son effect and also the Shapiro steps across the Josephson 
junction J J. For the ac Josephson effect, we assume that 
the JJ is shunted by a vanishingly small parallel resistor R\ 
and is biased by a large total current I to t ^> Ic (where I c is 
the critical current of the JJ) so that the DC voltage across 
the J J is Vdc ~ Itot/Ri- The finite voltage Vdc across the J J 
creates an ac current through the JJ, which is mostly car- 
ried by the large capacitor C and does not affect the voltage 
across R± . This ac current is inductively coupled to the large 
resistor R2 through the transformer T. By choosing a fil- 
ter Zf with a resonance at a frequency cj, one can measure 
the power-spectrum of the current (Isc(w)Isc(—w)), which 
would have peaks dependent on Vdc- This is the ac Joseph- 
son effect. For the Shapiro step experiment, Ri is chosen to 
be larger but still in the overdamped regime. An ac voltage 
v ac (t) = v ac ,o cos out is then applied (after making the resistor 
R2 vanishingly small). The dc voltage Vdc when measured as 
a function of Itot shows steps or kinks which are referred to 
as Shapiro steps. 

the voltage Vdc leads to a precession of the phase vari- 
able qb according to the equation with <j>(t) = AenVdct- 
This results in an ac current across the JJ 



Shapiro steps in terms of the radiated 
power-spectrum 

Introduction 

In this section we consider the Shapiro step experiment 
and compare it with the ac Josephson effect using the cir- 
cuit shown in Fig.[3j The ac Josephson effect occurs when 
a finite DC voltage Vdc is applied across the Josephson 
junction J J in Fig. [3j by setting the resistor Ri — >> and 
the current I tot correspondingly to 00. These two limits 
are equivalent to applying a nearly ideal voltage source 
V dc — hot/Ri- For the conventional Josephson junction 
with a current-phase relation 

J = / c sin0, (22) 



Ij{t) =I C sin (2eV dc t), (23) 

which is inductively coupled to the resistor R2 through 
the transformer T. The power dissipated in the resistor 
R2 as a function of the resonance frequency u of the filter 
Zf shows a peak at u = 2eV^ c , which is referred to as 
the ac Josephson effect. 

The Shapiro step experiment is in a sense the reverse 
of the ac Josephson experiment. An ac current is now 
injected into the transformer T in Fig. [3] through the 
voltage source v ac = v aCj o cos (jut) and the DC voltage V^ c 
is measured. To avoid externally imposing the voltage 
V^c ? one must now make the resistor Ri finite (in contrast 
to being vanishingly small) but small enough so that JJ 
is in the over-damped regime. To avoid extra dissipation 
R2 is made vanishingly small. The shunt capacitor C is 
now important to avoid having extra ac voltage arising 
from the ac current in the J J flowing through the resistor 
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R\. The voltage V dc is now given by the equation 

V dc = (I to t-Isc)Ri. (24) 

In addition, the transformer (assuming it to be ideal and 
1 : 1) induces an ac voltage — v aCj o cos ujt across the JJ. 
The minus sign is a result of the fact that the trans- 
former is in series with the junction and together these 
are shunted by a capacitor C, which forces the total ac 
voltage to be zero. The total current through the JJ is 
given by 



Isc(t) = I c sm{2eV dc t 



Vac,0 



UJ 



i I c sin (2eV dc t) - J cos (2eV dc + uj)t 
Zuj 



(25) 



The last term in Isc(t) is the only low-frequency term, 
which can lead to a low-frequency ac component Sv(t) in 
the voltage across the capacitor of the form 

Sv(t) = -I c Vac >° Rl — 

2uj x /1 + (2eV dc - uj) 2 RjC 2 

cos {(2eV dc - u)t + 4>(u)} (26) 

where (f)(uj) = tan -1 (R 1 C(2eV dc — lj)). This leads to an 
additional phase-fluctuation 

S<t>{t) = -IcRi / QC '° — — 

2uj(2eV dc - uj)^l + (2eV dc -uj) 2 R 2 C 2 

sin{(2eVd c -cj)t + 0H}. (27) 

Because of the pole when V dc = uj/2e, the above expres- 
sion is only small and meaningful as a perturbation for 
\2eV dc — lj\ > I c Ri- In this case, one can expand the 
low- frequency part of the expression for Isc(t) as 

SIsc(t) -tc^f cos {(2eV dc - uj)t + 5<j>{t)} 

~ ^-^f COS {( 2eV dc ~ Oj)t} 



I c Ri 



u ac,0 



^ 4uj 2 (2eV dc - oj) y/l + (2eV dc - uj) 2 R\C 2 
sin {(2eV dc - uj)t + <p(u)} sin {2(2eV dc - uj)t} 



(28) 



I c ^j- cos {(2eV dc - uj)t} 



J ac,0 



4to 2 (2eV dc -to) 



;{2(2eV dc -cv)t + <p(oj)} 



ilRi- 



u ac,0 



, (29) 



■ 4uj 2 (2eV dc - lj) y/1 + (2eV dc - uj) 2 R\C 2 ' 
The dc voltage in this case is determined by the equation 



hotRi — V dc + 



J ac,0 



{icRif 



4w 2 (2eV dc - uj) v / l + (2eV dc -ujyR 2 1 C 2 ' 

(30) 



The above equation has no solution in the range 



\T E> W K V ac$(I C Rl) 

\hotni- — \ > 



^ V 2e^l + (2eV dc -uj) 2 R 2 C 2 ' 

(31) 

which corresponds to voltage V dc very close to a;/2e, 
where the perturbation theory is expected to break down. 

In the range where V dc « uj/2e the voltage readjusts 
the phase, so that V dc adjusts to V dc = uj. This is the 
state which produces the Shapiro step in the current 
Itot\3Q\- In this case, the current in Eq. 25 has a DC 
component given by 



Isc,dc(t) = i c ^^COS( 



'rel 



(*). 



(32) 



where (j) re i(t) is the slowly varying (or constant) rela- 
tive phase between the superconducting phase and the 
ac perturbation. If (j) re i(t) remains bounded in fluctua- 
tions i.e. 4>rei(t) <C 27r, then it is referred to as the locked 
state. Since Isc, d c is bounded by I c ^f~ in magnitude 
the locked state can only be stable in the range 



UJ 



2eR x 



lr.— < Hot < 



Vac,0 



2uj 



2eR x c 2uj ' 



(33) 



where V dc plateaus at the value uj. The relative phase 
(t) is now locked to 



COS( 



(*) 



2eRx 



(34) 



In this range the voltage V dc shows a plateau in the neigh- 
borhood I to t ~ with width proportional to v ac ^. This 
range is approximately (i.e. overlaps only slightly) com- 



plimentary to the range in Eq. 31 where V dc has an os- 
cillating in time component. 



Classical fractional Shapiro effect 

Let us now consider how a doubled voltage Shapiro 
step can exist in the classical SNS junction with an 
avoided crossing. Using the classical Markov model in 
Eq. [7] ,we consider a current phase relation 



Isc(t) = s(f)I c sin(0(t)/2), 



(35) 



where s(t) = ±1 is a state variable that randomly fluc- 
tuates on a characteristic time T. If s(t) is constant over 
all time, based on the discussion in the previous section 
we expect that V dc is locked to V dc = uj/eR\ for a range 
of currents I tot ~ ujjeR\ which is given by Eq. 33 The 
DC current in the junction is now given by 



LSCAc 



Vac,0 
1 2UJ 



— S(t) COS <j) re i(t). 



(36) 



When s(t) changes from s(t) = +1 to s(t) = —1, the 
DC current in Eq. [36] flips sign. This results in the DC 
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voltage jumping to ^ - 2I c R ± ^f cos <j> re i(t) . The DC 
voltage is now off resonance and leads to a change of 
(t) ~ — 2I c Ri^r^- cos (j) re i(t). Based on this equation, 



»rel\ 



one expects the system to fall back in lock on a time-scale 



V ac,0 



2uj 
I C R\ 



(37) 



If the T i.e. the time between flips of s(t) satisfies the 
constraint T ^> r, then one expects plateaus at Vdc = w/e 
corresponding to the doubled voltage Shapiro steps. On 
the other hand, the limit of small v ac is behavior is some- 
what different and one expects only kinks in the Vdc ver- 
sus I tot curves instead of true plateaus. This is because in 
the weak v aCy o limit, we expect the voltage V dc at a fixed 
current bias I to t to be analytic in v aCy o and also indepen- 
dent of the phase reference of the microwave. Therefore 
we expect Vdc to vary with v ac o with the second power 
v 2 

The presence of the step in the voltage as a function 
of current Vd c {hot) leads to a divergence in the function 



f(V dc = ItotRl) 



hotRl — Vdc 

Ri 



(38) 



The slight subtlely in the abve definition, which is that 
ItotiYdc) is n °t a unique function of Vd c (Itot) is technical 
and is avoided whenever the plateau has a slight slope. In 
the small v aCj o limit we expect all quantities, including 
the function / to be finite (corresponding to I to t as a 
function of Vd c being smooth) and proportional to / ex 
v 2 c . The fact that at larger values of v aCy o a plateau in 
Vdc(hot) (implying a divergent /) suggests that f(V dc ) has 
a kink at the values of voltage where there is a plateau 
in the noiseless J J limit. We will discuss how such a 
kink in the function / is related to the peaks in the noise 
spectrum in the v aCj o —> limit in the next section. 



Shapiro kinks in the noisy Josephson junction 

The small v aCi o limit in a noisy Josephson junction of 
the type considered in the previous section, in general 
does not give rise to a linear in v aCy o dc voltage plateau 
or a full phase lock. Instead it gives rise to phase cor- 
relations that result in a kink with strength v 2 c as we 
discuss in this section. We will find that the Shapiro 
kink is a more universal feature of the Josephson effect 
and depends on fewer of the details of locking. 

Let us consider the circuit in Fig. [3] in the Shapiro step 
regime so that we can set R2 = and R\ is finite. Let 
us define the integral 



dt^Vdcih) - ItotRi) 



(39) 



the field <j>(t) and hence Vdc across J J can be calculated 
through a Keldysh action [34] 



S[(j) ch (j) q ] = 
-2 J dt</> q \ 



Cd 2 <\> cl , 



1 d(j) c i 



+ i 



2R ± 



2i 



dtdt' 



dt 2 2R X dt 
2smh 2 (irT(t-t , ))/7rT 2 ' 



(40) 



where (j) c i and <\> q are the classical and quantum parts of 
the fluctuating field 0, T is the temperature of the resis- 
tor, 5q is the action of the J J after the ABSs have been 
integrated out. We will now consider the high temper- 
ature limit where the quantum (or shot) noise (i.e. the 
last term) can be neglected compared to the Johnson- 
Nyquist noise and moreover, we will assume that T is 
large enough so that <p q is small and we can expand the 
action to lowest order in 6 n . 



S[<f>ch<l>q 

T 



dt<j) t 



r C d 2 (f) c i ^ 1 d<j> c i 
2R 1 dt 



(i Wi +d 2 q S [^ci:0])^ q . 



2 dt 2 

2 



(41) 



Following Ref. [51] . the correlators of the classical field 
<pci (t) can be obtained by integrating out <p q and replacing 
it by a Langevin equation for <f> c i(t) 



C 



1 



2 dt 2 

Ri 



2Ri dt 



iscttdKt) = at) 



2T 



dl q s [4>ci,o](t,t'nt') 



(42) 



where is Gaussian correlated noise, which in the 
Nyquist noise dominated (i.e. high temperature) limit 
has a correlator 



(H(tW)) = 5(t-t>) 



T 
R~i 



(43) 



To obtain Eq. 42 we have expanded to lowest order in 
d 2 So, which is the noise spectrum of the current in the 
J J i.e. we have made the assumption that the noise from 
J J produces only a small voltage across R compared to 
V dc - The function Isc[<l>ci](t) = <9</> q So[0d, 0q]\<f> q =o is the 
expectation value of the supercurrent in the presence of 
0czWj which is expected to vanish in the dc limit. Eq. 42 
is written more explicitly in terms of the noise correlator 
as 



Cd 2 ^ 



dt 2 
Ri 
2T 



d t 1 d(j> c i 
2R 1 dt 



f 

J — ( 



dt'{I(t)I(t')) c [cl> cl (ti)W), 



(44) 



as the phase fluctuation variable associated with the volt- 
age difference across the capacitor. The correlators of 



where (I(t)I(t')) c [(j) c i(ti)] is the connected part of the 
current-fluctuation correlator from J J. 
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Defining (t) to be the solution to Eq. 42 for S = 



i.e. 



A (0)'" ' -!*=^ 



'(*) 



-e-^^)^(i'), (45) 



Eq. 44 may be rewritten as 



/ ^! e -^ /" df</(t")j(f)> c [0c J (ti)]^(f). 

J — oo J —oo 

(46) 



The mean dc voltage shift as a result of the junction J J 
is given by 



(7(t")f(t')M#S"(*.)]«*') 



(t-t") 

RC 



t" 



dt' 



(47) 



In the small i? limit the variation of the phase is slow so 
that we can write 

M Ri t l dt" ct-t") /•*" , , 

/g^™ w ,, (ti) _ ,, (tW)i (48) 

J ^ (*l) 



Substituting Eq. [45] and averaging over the Langevin 
noise leads to 

*Vdc(t) = ^ / dt"dt'dt x e-^r(\ - e"^H 

Jf <t x <t" <t 



I 



2C . 

d<f> S(I(t")I(t')) c 
2- i^) (tl ) 



(49) 



The integral over the phase <\> = </>^(£'), which represents 
the part of the phase that diffuses because of the thermal 
noise from averages over all initial phases. The shift 
in DC voltage 5Vd c — Vd c — ItotRi is directly proportional 
to the Shapiro kink function / defined in Eq. |20j 

Eq. [49j which relates the shift of the voltage associated 
with a Shapiro kink to a derivative of the current-current 
correlator is the result of the derivation in this section. 
Because of the integrals over time the DC voltage Vd c , 
contains only contributions from the low-frequency part 
of I(t) in the current-current correlator, which are pro- 
portional to v ac ^. Therefore SVd c (t) v 2 c as conjec- 
tured. Furthermore, the size of this response provides the 
size of the kink (apart from the factor of v 2 ac ). Moreover 
we expects peaks in the noise spectrum, which is what is 
measured in the ac Josephson effect to correlate with the 
steps in the response of the noise. 



Computing the Shapiro kink 

The size of the Shapiro kinks are proportional to the 
derivative of the connected correlator (I(t)I(t')) c , which 
must also be expanded to second order in v ac . Therefore 



5{I{t")I{t')) . 



is a 



the relevant low- frequency component of 

third order derivative of (I(t)I(t')) c with respect to the 
phase. Considering the expression Eq.[l7j we notice that 
terms where the derivative with respect to phase <p op- 
erates on the U terms are suppressed by factors propor- 
tional to the critical current associated with the ABSs. 
Therefore, we can expect the dominant contribution to 
the derivative to come from the term where all of the 
derivatives with respect to <p act on the phase-dependence 
of the current operator itself i.e. 



S(I(t")I(t')} ^ vl 



6(t" - h) 



,.S 2 I(t")SI(t'), 



[COS 2 Ut" ( " 'JC ' I (t' ) ) c + COS Ut " COS Ojt ' ( - '- ) c 



■ cos cut { 



„<*/(*") 5 2 I(t') 



)c] + A, 



(50) 



where A contains the rest of the terms which have higher 
order correlators of the current /. 

The second order in current / contribution to Vd c is 
thus written as 



6Vdc W ~ 2Cw 2 



r (t- t ") 
/ dt"dt'e~^^(\ 
Jt'<t"<t 



/ 



^<™ 2 urf\l"\i!')l(i!)) c + cM 



2tt 



+ cos" ut'{I'{t")I"{t')) c ), 



(51) 



where J'(0) = d+I (</>), = d 2 I($) and J"'(0) = 

<93/(0). Rewriting the products of the cosines as sums 
of cosines, we observe that since both t' and t" are 
integrated out, only the second term can have a low- 
frequency contribution that is significant, so that 



&v£> i')^'.'':; / dfdt'e- — ( 1 - e" — ) 

UJ Jt'<t"<t 



COSCJ 



(t"-f) J ^(I"(t")I'(t')) c 



(52) 



where the average over the overall phase <p is implicit in 
the remaining equations. 

Following Eq. [lT| the connected correlator for a pair 
of operators 0^\t) and 0^ 2 \t) is written as 



(O (1) (t")O (2) (0)<: ~~ ^ n{1) ^" v ° (2) * (+f 



Ai ,A2 



-T\t'-t"\ i( €Xl -ex 2 )(t'-t") 



(53) 



where 0^ Xe (t) = (ip Xb {t)\0^\t)\^ Xc (*)> are r-periodic 
functions. Assuming T > 1/(R\C) i.e. the broadening of 
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the peaks appears is intrinsic, 

sv. 



1 ir f (*-*") 

d u J dt'smujt' J ^(I"(t")I'(t' + t")) c . (54) 

Substituting the expression for the connected correlator 
from Eq. [53l leads to the result 



0V dc 4C UJ 2 ^ i Ai,A 2 W i Ai,A 2 l 

[ 1 

l {uj-iT + (e Xl -e X2 +nftj)} 2 

1 

+ {u + ir-(e Xl -e X2 +nSlj)y 1, 



n) 



(55) 



where I' XiXa (n) = £ dtj na < t (<p Xl W'(t)\<Px a (t)) and 

'W") = Iodte in ^{ip Xl (t)\I"(t)Mt))- Note that 
averaging over </> has no consequences in the high T and 
low frequency limit used for this calculation. 



The above dc voltage has peaks at the bias voltages 
I tot, where the denominators e\ 1 — e\ 2 + nVtj — uj vanish 
can be expected to show a kink structure. However, the 
induced dc voltage 5Vd c has a profile which is related to 
the derivative of the peak-structure corresponding to the 
ac Josephson effect shown in Fig. [2] 



